Formation of subgap states in carbon nanotubes due to a local transverse electric field 



Jesse M. KindeiQ and E.J. Mele 
Department of Physics, University of Pennsylvania 
Philadelphia, PA 19104 
(Dated: February 2, 2008) 

We introduce two simple models to study the effect of a spatially localized transverse electric 
field on the low-energy electronic structure of semiconducting carbon nanotubes. Starting from the 
Dirac Hamiltonian for the low energy states of a carbon nanotube, we use scattering theory to show 
that an arbitrarily weak field leads to the formation of localized electronic states inside the free 
nanotube band gap. We study the binding energy of these subgap states as a function of the range 
, and strength of the electrostatic potential. When the range of the potential is held constant and 

the strength is varied, the binding energy shows crossover behavior: the states lie close to the free 
nanotube band edge until the potential exceeds a threshold value, after which the binding energy 
q , increases rapidly. When the potential strength is held constant and the range is varied, we find 

resonant behavior: the binding energy passes through a maximum as the range of the potential is 
increased. Large electric fields confined to a small region of the nanotube are required to create 
\Q , localized states far from the band edge. 
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PACS numbers: 73.20.At, 73.22.-f, 73.63.Fg, 73.63.-b 



Single- walled carbon nanotubes are cylindrical fullerenes with a radius on the order of a nanometer and lengths that 
can exceed a micrometer. They occur in both semiconducting and metallic species. The conductivity has a purely 
geometric origin that allows all nanotubes to be divided into three types: moderate-gap semiconductors, narrow gap 
semiconductors, and metals. The most common is the first type, which have a band gap between 0.1 eV and 1 eV 
that scales inversely with the tube radius R. This band gap is required in order for the electron wave function to be 

S single- valued. The narrow-gap semiconductors have a band gap on the order of 0.01 eV that scales inversely with R 2 
and arises from the distortion of the honeycomb lattice when it is wrapped into a cylinder. In the metallic armchair 
i-^j , nanotubes, the formation of a band gap is prohibited by symmetry^ We only consider moderate-gap semiconducting 
c"j ' nanotubes in this paper. 

In applications involving carbon nanotubes as electronic and optical components, a method for controllably and 
continuously varying the size of the intrinsic band gap would be useful. In principle, elastic deformations and magnetic 
flux through the a nanotube can be used to control the band gapJu£ Another possibility is to tune the band gap with 
t-H ■ a static electric field perpendicular to the nanotube axis. A transverse electric field leads to a Stark shift of the 
electronic energy levels. The scale of this shift is set by the electric potential difference across the diameter of the 
nanotube: AV ~ 2EqR, where Eq is the magnitude of the electric field. In analogy with the Stark effect in atomic 
physics, one might expect a continuous reduction in the nanotube band gap as the strength of the electric field is 
increased. This is the case for boron-nitride nanotubes, which are wide-gap semiconductors with a band gap on the 
order of 5 eV. A Stark shift was predicted and has been observedi^^ 

Several authors have studied the effect of a uniform static transverse electric field on the band structure of moderate- 
gap single-walled carbon nanotubes 1 1 1 1 1 1 2 Instead of a linear or quadratic Stark shift in the band gap, nu- 
merical simulations indicate crossover behavior. The change in band gap is small until the electric field exceeds a 
threshold, after which the band gap decreases linearly with increasing field strengthi^^ ' 10 ! 11 

Novikov and Levitov have used the low-energy Dirac theory for carbon nanotubes to provided a natural explanation 
of these numerical results^ In the Dirac theory, the band gap is fixed at its zero-field value until the electric field 
exceeds a critical strength. Below the critical field, the band gap is protected by a chiral gauge symmetry; above the 
& \ critical field, the effective mass of the charge carriers changes sign, and the band gap closes around the symmetry- 
protected point. 

Numerical simulations do not exhibit critical behavior — there are small increases or decreases in the band gap, 
even in very weak fields. However, changes in the band gap are small below a threshold field strength. As shown in 
App.El the critical field predicted by Novikov and Levitov is smaller than the threshold fields reported in tight binding 
calculations. The work of Novikov and Levitov shows that the linear Dirac theory provides a qualitative description of 
the response of semicondcuting nanotubes to a uniform transverse electric field: crossover behavior. This simple model 
captures the essential physical properties of a nanotube necessary to understand the general features its electronic 
response. 

In this paper, we study a nanotube in a nonuniform field and explore the possibility that spatially varying fields 
could couple to the minimum gap. The chiral symmetry of the uniform field model is broken by an inhomogeneous 
electric field, and the electronic response to a nonuniform field is qualitatively different. 
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Our work suggests the formation of exponentially localized subgap states is a general feature of the electronic 
response of semiconducting carbon nanotubes to transverse electric fields localized to a small region along the tube 
axis. Since the states are localized, they do not reduce the band gap for free particle transport along the nanotube 
axis. This is different than the response to a uniform field, in which the gap for excitation from the valence band to 
the conduction band is modified. Our results indicate that this gap remains fixed in a field localized to a small region 
along the nanotube axis, even when the field strength is large. 

We study two simple models for a localized electric field. In both models, the potential modifies the density of 
states in two ways: 1) electron scattering leads to periodic variations in the density of states along the nanotube axis, 
and 2) the potential creates exponentially localized states inside the free nanotube band gap. Fig. [I] is a plot of the 
density of states near the band edge for a free nanotube and the two models we studied. The potentials define the 
center of an infinite nanotube. The plot shows the density of states at a distance of one tube radius from this origin. 




FIG. 1: (Color Online) The local density of states near the center of a carbon nanotube for energies close to the free nanotube 
band edge at e — 1/3. Energies are given in units of hvp/R. The plot shows bound states inside the band gap and a modification 
of the density of states near the band edge. 



Localized states form inside the free nanotube band gap in arbitrarily weak fields; however, the binding energy as a 
function of potential strength shows crossover behavior similar to that observed in numerical simulations of nanotubes 
in uniform transverse electric fields. When the range of the potential is held constant, the subgap states lie close to 
the band edge until the electric fields exceeds a threshold strength. 

The remainder of the paper is divided into four sections. In Sec. HI we present the Hamiltonian for a semiconducting 
carbon nanotube in a static transverse electric field and outline the scattering theory we used to study the spectra 
of these systems. Next, we analyze two models models in which the potential due to the applied field is localized 
along the nanotube axis. In Sec. |Tl] we consider an electric field that can be modeled as a delta function in the 
long-wavelength Hamiltonian. In Sec. IIIIl we consider an effective square well potential. In Sec. II V[ we compare the 
solutions of these models. Our main results are presented in Figs. [3l |H and [5] 



I. HAMILTONIAN AND SCATTERING THEORY 



The electronic structure of a carbon nanotube is inherited from the band structure of two-dimensional graphene, 
a planar honeycomb lattice of carbon atoms. Graphene is a semimetal with a vanishing band gap at the six corners 
of the first Brillouin zone, the K and K' points^ Since the honeycomb lattice has a two atom basis, the Bloch wave 
functions are a sum of amplitudes for each sublattice, which we denote A and B: 

#k(r) =Vu,k(r)+Vs,k(r). (1) 

Near the K and K' points, the electron wave function can be approximated by introducing envelope functions 
to describe long-wavelength variations of the exact wave functions at K or K'. Thus, if the total wave vector is 
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k = K + q, then 

*k(r) = <Pa,k( t ) "kW + ^s,K(r) w K {r). (2) 

The Bloch functions ip a are the exact K wave functions and oscillate on the scale of the lattice spacing. The envelope 
functions v and w oscillate on much longer length scales. 

When expanded to linear order in q, the resulting effective Hamiltonian for the envelope functions near the K 
point 19 is a massless Dirac equation? 1 - ' 14 ! 15 

- *** + few) = £ w ■ (3) 

The factor /w^ comes from the tight-binding model of graphene. Its value is 3at/2 0.53 cV • nm where t is the 
nearest neighbor hopping energy and a is the distance between neighboring atoms. The effective Hamiltonian describes 
spinless electrons in the graphene lattice. The two components of the pseudospinor give the relative amplitude for an 
electron or hole to be localized on the A and B sublattices. 

The low-energy theory of a carbon nanotube follows from this Hamiltonian. One defines a circumference vector C 
in the graphene lattice, which determines both the radius R and chiral angle 9 of the nanotube. The full electron wave 
function must be single- valued on the surface of the nanotube: 4'k(r + C) = x f'k(r)- This requirement collapses the 
energy surface of graphene into a set of discrete energy bands. If none of the energy bands passes through the K and 
K' points, then the nanotube will have a band gap inversely proportional to its radius, and the envelope functions v 
and w satisfy quasiperiodic boundary conditions.' 1 " These are the moderate gap semiconducting carbon nanotubes we 
consider in this paper. 

The natural description of a nanotube is in cylindrical coordinates. The displacement along the nanotube axis, z, 
and angular displacement around the circumference, 4>, are related to the (x, y) coordinates of the graphene sheet by 
a rotation through the chiral angle 9. (We define 9 to be the angle between the K point and the nanotube axis, so 
9 = for armchair nanotubes and 9 — tt/6 for zig-zag nanotubes.) The transformations are 



z \ _ / x cos 9 + y sin 9 
R4>) \ycos9 — xsm.9 



(4) 



In the coordinates of the nanotube, Eq. {SJ can be expressed 

H = -i{hv F /R) e - 1 ^ 9 ' 2 (Ra x d z + Uyd^) e ia " e l 2 . (5) 

The energy spectrum of this free nanotube Hamiltonian is a series of hyperbolic bands given by 

£ m (q) = ±(Hv F /RW(qR) 2 +Al, (6) 

where q is momentum along the nanotube axis and A m = m+ v/3. The integer m is the subband index, and v = ±1 
is required by the quasi-periodic boundary conditions on the envelope functions. Without loss of generality, we choose 
v = 1 in the remainder of the paper. 

Placing the nanotube in a static electric field perpendicular to its axis introduces the scalar potential 

V(z^)=e(E /K)Rf(z)co S( f>. (7) 

Eq sets the scale of the applied field strength; variations in field strength along the nanotube axis are described by 
the dimensionless function f(z). The electric field that enters the Hamiltonian is the screened electric field, not the 
applied electric field. Calculations indicate the effects of depolarization are significant and reduce the strength of the 
external field by a factor of roughly 5.— ^ We include depolarization effects through the constant n. 

The natural energy scale of the nanotube system is Hvf/R- We will use this quantity as the unit of energy and 
work only with the reduced Hamiltonian. The reduced potential strength is u = eE$R 2 / nhvp. This is a convenient 
parameter to use in calculations, but it obscures the magnitude of the fields involved. For reference, u = I in a 
nanotube with a radius of 1 nm corresponds to an applied electric field strength of about 10 V/nm. 

The reduced Hamiltonian — with dimensionless eigenvalues e — is the starting point for both of our models: 

H = -ie- la * e/2 {Ra x d z + ayd^) e la * e l 2 + u f(z) cos </>. (8) 



Because the potential depends on the azimuthal angle </>, angular momentum is no longer a good quantum number. 
However, the potential is single-valued around the circumference of the nanotube. The Hamiltonian has a discrete 
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symmetry under <j> — > <f> + 2ir, so its eigenstates can be expanded in the basis of angular momentum eigenstates of 
the unperturbed Hamiltonian. The potential is proportional to cos </>, so it will only connect states whose angular 
momenta differ by ±1: i.e., electronic states in neighboring energy bands. Projecting out the azimuthal degree of 
freedom in the Hamiltonian of Eq. © gives 

H m , n = e-^ 6 ' 2 (a y A m - iRa x d z ) e^ e / 2 5 m , n + (u/2) f(z) 6 m , n±1 , (9) 

where m and n are the band indices of the unperturbed nanotube. 

We study this Hamiltonian for two different models of a nonuniform transverse electric field. In the first, we treat 
the potential as a point scatterer with /(z) oc S(z); in the second, we model the potential as a square well so that f(z) 
is constant within a region of length L and zero outside. Neither of these models can represent a real electric field 
because neither potential satisfies Maxwell's equations. However, common features in both these models suggests a 
universal response to a localized static transverse electric fields: the formation of localized states with energies inside 
the unperturbed band gap. 

Eq. ((9J is an effective Hamiltonian for the long-wavelength modes of the carbon nanotube. Accordingly, f(z) is an 
effective profile of the electric field. The physical interpretation of f(z) can be made precise by considering the three 
important length scales in the problem: the length over which the potential changes significantly, Ax; the wavelength 
(or decay length) of the envelope functions, A = 27r/|q|; and the distance between neighboring atoms of the honeycomb 
lattice, a. The natural cutoff for the linear Dirac theory of electrons in carbon nanotubes is A ~ R. 20 Particles with 
energies near the band edge may have wavelengths much larger than the radius, on the order of the tube length. As 
a result, we may consider the wavelength of the envelope functions to be bounded from below: A > R. 

We consider potentials that satisfy a <C Ax <C A. The physical field is continuous. The change in potential over 
distances comparable to the lattice spacing is negligible, but on length scales comparable to the wavelength of the 
envelope functions, the field strength changes abruptly. Since the physical potential is roughly constant over distances 
comparable to the lattice spacing, there is no intervalley scattering between the K and K' points even though the 
potential in the effective Hamiltonian is discontinuous. For such a potential, the energy spectrum can be determined 
by considering a single valley in isolation. Fig. [5] illustrates the types of potentials to which the simple models discussed 
in this paper might apply. 

In the remainder of this paper, we will ignore the chiral angle 9. The dependence of TL m .n on 9 can be removed 
with a unitary transformation. This is discussed in App. [Bl in which we show that the density of states and bound 
state energies of Eq. ^ do not depend on the chiral angle. The results derived for 8 = apply to chiral tubes as 
well.21 




(a)Point Potential (b)Square Well 



FIG. 2: (Color Online) Field profiles that can be analyzed with the models presented in this paper. The black lines show the 
electric field strength along the nanotube axis. The red lines indicate spatial variations in the envelope functions along the 
axis. The potentials are roughly constant on length scales comparable to the interatomic spacing, but change significantly on 
the scale of the envelope function wavelength. 

In (a), the potential rises from zero to its maximum value and falls back to zero within a single wavelength of the envelope 
function. In the effective Hamiltonian, this potential can be modeled as a point scatterer. In (b), the potential rises from zero 
to its maximum value over a distance smaller than the envelope function wavelength. It is constant for several wavelengths, 
then falls back to zero over a distance smaller than the wavelength. In the effective Hamiltonian, it can be represented by a 
square well. 

To study the effects of these potentials on the electronic structure of a carbon nanotube, we use scattering theory to 
determine the Green function of Eq. §§§ from the Green functions of a free nanotube. The singularities of the Green 
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function give the energy spectrum, and its trace gives the local density of states. Free nanotube Green functions are 
shown in an open font: G. The full Green functions for a nanotube in the applied field are displayed in a script font: 
6. 

The Green function for a charge carrier in band m of a free nanotube that satisfies outgoing boundary conditions 
as \z — z'\ approaches ±00 is 

G ( Zj zl) = s + a y A m +s g n( Z - z')a x U(e) ^ lz . z y R 



where f m (e) = ^e 2 - A^. This expression is derived in App. [Cl The Green functions are 2x2 matrices to account 
for the psuedospin of the envelope functions^ When \e\ < A m , £ m (e) is imaginary and the Green functions decay 
exponentially as \z — z'\ increases. 

The Green functions of the perturbed system are solutions of Dyson's equation: © = G + <&V &. Formally, the 
solution is 

6 = — - — G, (11) 
where 1 is the 2x2 identity matrix. Expanding this expression gives an infinite series: 

e = G + Gvlj2 \ <G - ( 12 ) 

lra=0 J 

This is an operator equation, and products are shorthand for integration over position along the nanotube axis and 
summation over all bands. In general, this equation cannot be solved exactly. We make a series of approximations 
that allows us to carry out a partial summation and approximate the exact Green function. 

As shown in App. [D] when the energy of a charge carrier is close to the band edge, the largest terms in the Dyson 
series are those that only contain transitions to the neighboring bands. Retaining only Green functions from the first 
three bands (m = 0,±1), Eq. (fT2|) can be summed to give the full Green function for the lowest band: 

6 = Go + G r - * - G . (13) 

JL — HjoI 

To = V(G+i + G~i)V is an effective potential for the lowest band. It will be nonlocal in general; however, for the 
models studied in Sec. HT1 and IIIH the effective potential only depends on energy. The general form of the full Green 
function is then 

6 (e, z, z') = G (e, z, z') + G (e, z, 0) r (e) - - 1 . G (e, 0, z'). (14) 

1-Go(e,0,0) r (e) + i?/ 

A small imaginary number, irj, has been included in the denominator of Eq. (|14[) so the poles in the complex energy 
plane correspond to physical particles. This is the retarded Green function. 

The density of states in the nanotube is proportional to the trace of the Green function. The poles of 60 are bound 
state energies, and branch cuts correspond to energy bands i 17 ' 18 When the energy approaches a band edge of the 
unperturbed Hamiltonian, either Go or To diverges. In either case, 60 remains finite, which means 60 does not share 
any singularities of the free nanotube Green functions. (This is clear in Fig. [TJ The free nanotube density of states 
has a van Hove singularity at the band edge, while the density of states vanishes at the band edge in the presence of 
a potential.) 

If there are any poles in Eq. (|14[) . they occur when the determinant of the denominator vanishes: i.e., 

det[l-G (£,0,0)r (e)] = 0. (15) 

This equation yields a sixth-order polynomial, and we were unable to find a general solution for energy as a function 
of potential strength. To find the roots, we created a simple numerical routine to determine the energy (within a 
user-specified tolerance) at which the determinant changes sign. 

In both of the models, numerical results suggest that two localized states are created inside the band gap for any 
nonzero potential — one for the valence band and one for the conduction band. In App. [El we show analytically that 
the delta function model has states inside the free nanotube band gap for arbitrarily weak fields. The analog of this 
proof for the square well model is discussed in Sec. IIII1 

The fact that bound states come in pairs reflects the symmetry of the particle and hole spectra for the potential in 
Eq. 0. When the charge of the particle is reversed, V(z, <fi) — > —V(z,(j>), which is equivalent to physically rotating 
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the nanotube by 7r about its axis. Since the point at which <f> — is arbitrary, the particle and hole spectra must be 
identical, and all energy eigenvalues of Eq. © — including bound state energies — come in pairs ±e. 

In the more familiar case of a charged impurity, an attractive potential for electrons is repulsive for holes, and the 
bound state spectrum is not symmetric. The difference between the spectra of a charged impurity and the potential 
considered in this paper arises from the angular dependence of the latter. 

We define sb to be the (non-negative) energy at which the full Green function has a pole. To study the dependence 
of the bound state energies on the range and strength of the potential, we plot the binding energy 7 = 1 — e^/Ao, 
which is the normalized difference in energy between the localized state and the free nanotube band edge. Our 
numerical results show < 7 < 1. If 7 = 0, there are no states inside the gap; if 7 = 1, the subgap states lie on top 
of each other at the center of the free nanotube band gap. 

In the following two sections, we apply the general formalism described here to two specific models. In the first, we 
treat the effective field profile as a delta function. In the second, we model the effective field profile as a square well. 
We study the binding energy of the subgap states as a function of the potential size L and strength u. Figs.[3j[4l and 
[5] illustrate our results. 

II. DELTA FUNCTION MODEL 

In this section, we study the electronic structure of a carbon nanotube in a transverse electric field that can be 
treated as a delta function in the effective Hamiltonian. (See Fig. [5{a).) 

The physical potential can be written V(z) = Voh(z), where h(z) is some function that is localized to a small 
region of the nanotube axis — for instance, a Lorentzian distribution peaked around z = 0. We define the range of 
the potential to be L — J dzh(z), then replace f(z) by L8(z) in Eq. ©. This ensures that the integrated strength of 
the potential is unchanged. A delta function potential in the effective Hamiltonian will be a reasonable approximation 
as long as L < R. 

The effective Hamiltonian for the delta function model is 

W m ,„ = (<T y A m - iRa x d z ) 8 m , n + (uL/2) 5{z) S m ^ ±1 . (16) 

Since the potential is proportional to uL, the range L and strength u of the potential cannot be varied independently. 
Any particular choice of (u, L) is equivalent to an infinite class of choices (it/3, Lj ' (3). 

When the potential is proportional to a delta function, the integrals of Eq. (fT2|) can be evaluated exactly. For 
instance, the term < GjqV < &±iV < &o is 

{uL/2) 2 G (e, z, 0)G±i(e, 0, 0)G (e, 0, z'). (17) 

As a result, the full Green function is given by Eq. (frjj with 

r (e) = (uL/2) 2 [G +1 (e, 0, 0) + G_i(e, 0, 0)]. (18) 

The local density of states is proportional to the imaginary part of the trace of Eq. (fl4|) . When e lies inside the 
free nanotube band gap, the Green functions are exponentially decaying functions, as mentioned in Sec. fl] Since the 
bound states lie inside the free nanotube band gap, the local density of states decays exponentially with increasing 
distance from the potential: 

P (e B ,z)cxexp[(2\z\/R)yjA 2 -e 2 B ]. (19) 

The bound state wave functions are centered on the potential and decay exponentially, like the solutions of the ID 
Schrodinger equation with an attractive delta function potential. 

Fig. [3] shows 7 = 1 — £s/Ao as a function of potential strength u when L = R. The binding energy shows crossover 
behavior similar to the band gap in a uniform transverse electric field. Although the subgap states are very near the 
band edge when the field is weak, there is no critical behavior like that found by Novikov and Levitov. In App. [El 
we show analytically that the Green function in Eq. (|14[) has two poles for an arbitrarily weak potential and that the 
binding energy of the subgap states increases as u for small u. 

This quartic scaling is not unique to the Dirac model and occurs in nonrelativistic models as well. If the potential of 
Eq. ([7]) is introduced to a Schrodinger equation describing electrons confined to the surface of a cylinder, the binding 
energies also grow with u 4 . Because the potential only connects states in neighboring bands, the effective potential 
in a particular band is proportional to u 2 8{z). The binding energy will scale with the square of the effective potential 
strength, which explains the u 4 dependence in Fig. [3J 

The figure shows that the potential must be strong (cSqL > hvp/R) to produce subgap states that lie more than 
a few percent below the band edge. 
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FIG. 3: (Color Online) Binding energy as a function of potential strength. The horizontal axis shows the value of the 
dimensionless parameter u — eEoL/(hvF / R) that characterizes the applied field strength. The line is 7 = 1 — es/Ao for L — R. 
The subgap states approach the middle of the gap as the potential strength is increased, but lie close to the band edge in weak 
fields. 



III. SQUARE WELL MODEL 

In this section, we model the potential as a square well in the effective Hamiltonian. (See Fig.[2jb).) In this model, 
the strength and size of the potential can be varied independently, and the restrictions on the range of the potential 
are relaxed. We find the effect of changing the range of the potential is different than changing the field strength. 

A square well potential in the effective Hamiltonian of Eq. @ corresponds to a physical field that increases from 
zero to its maximum strength in a distance that is small compared to the scale on which the envelope functions vary. 
The field strength is assumed to be constant for some length L, after which it decreases back to zero. 

For the square well model to be a good approximation of the physical potential, L must be smaller than the length 
of the nanotube. Otherwise, the electric field should be treated as uniform. There are no other physical restrictions 
on L. However, we are only able to solve Dyson's equation after imposing the additional restriction that qL <C 1 
where q is the wave vector of the lowest band (m = 0). Since we study states near the band edge where q approaches 
zero, this is not a restrictive assumption. It makes the effective potential a point scatterer for states in the lowest 
band. The difference between this model and the delta function model is that the (exponentially decaying) states in 
neighboring bands (m = ±1) interact with the potential in a region of finite size. 

The effective Hamiltonian for the square well model is 

H m ,n = (o-yAm. - iRa x d z ) 5 m ,n + (u/2) Q(z) Q(z - L)5 m ,n±l, (20) 

where Q(z) is the step function. 

In the delta function model, 8(z) set all of the nonanalytic functions of the Green functions in Eq. (fTU)) equal to 
zero in the integrals of the Dyson series. For the square well model, the integral has to be evaluated piecewise to 
account for the terms in the Green function proportional to sgn(z — z') and \z — z'\. For instance, the second-order 
term 

G VG ±1 VG = {u/2) 2 / dx dyGo{e,z,x)G ±1 (e,x,y)G {e,y,z') (21) 
Jo Jo 

must be divided into eight regions of integration if both z and z' lie inside the square well. We are unable to derive 
a general formula for the 0(u 2n ) terms of the Dyson series. As a result, we cannot sum the series, calculate the full 
Green function, find the density of states, or search for localized states below the band gap. 

By requiring that qL <C 1 where q is the wave vector of the lowest band, we can approximate the integrals and 
sum the series of Eq. p2[) . This approximation simplifies the expression for Go(e, z, z') in two ways. First, if L > R, 
then qR is smaller than qL, which is negligible. Thus, we can ignore the term in the Green function proportional 
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to sgn(z — z'). Second, when Go appears inside an integral, its phase, exp(iq\x — y\), can be replaced by a more 
convenient function. If x and y lie inside the region of integration and z is arbitrary, the error introduced by the 
following substitutions is of order qL: 

e iq\z-x\ 

e iq\x-y\ > e iq(x+v)_ (22) 

(Since z may lie anywhere along the length of the nanotube, we do not assume q\z\ is small.) As an example, consider 
the following: 

pL pL pL pL 

/ / dxdy e lqlx - yl = / dxdy e lq{x+y) +0{qL). (23) 
Jo Jo Jo Jo 

Since we are interested in the limit qL <C 1, the error is insignificant. 
Assuming qL <C 1, the integral in Eq. (f2Tj) gives 

G VG ±1 VG a = u 2 A±i(e) G (e, z, 0)G±i(e, 0, 0)G o (e, 0, z') + O(qL), (24) 

The scalar function A±i(e) depends on e and L. 

Aa(e) = 2{A f_ A§) (x a (e) + e-*^ - l) , (25) 

where x Q = (L/R)^A 2 a -e 2 = i&(e) [L/R). 

Ignoring the O(qL) terms, , all the integrations and summations can be carried out as in the delta function model 
and the full Green function is given by Eq. (TT4"|) with 

r (e) -u 2 [A+i(£)G +1 (£,0,0)+A_ 1 ( £ )G_(£,0,0)]. (26) 

Comparing this expression with Eq. (|18p shows that when qL <C 1, the square well model is equivalent to a delta 
function model with an energy-dependent range: l a (e) = 2y/A a (e). As a result, the arguments of App. [El also apply 
to the square well model. 

When L <C R and e ~ A , A(e) approaches L 2 /4, and the square well model reduces to the delta function model, 
for which there are bound states at any field strength. When L ^> R, the exponential term in Eq. (|25p can be ignored. 
An analysis similar to that of App. [E] shows that there are bound states for arbitrarily weak fields in this limit as 
well. Since A Q (e) is a smoothly increasing function of L, we interpolate between these two limits and conclude that 
a square well potential of any strength and range leads to the formation of localized states inside the free nanotube 
band gap. 

As in the delta function model, the bound state wave functions decay exponentially away from the potential , and 
the binding energy is proportional to u A for small u. 

Fig. 2] shows the binding energy of the subgap states as a function of electric field strength when the range of the 
potential is fixed (fixed L, variable u). The shapes of the individual curves are similar to that of the delta function 
model in Fig. [3] The range of the potential determines the asymptotic value of the binding energy as the field increases. 

Fig. [5] shows the binding energy of the subgap states as a function of potential width when the potential strength 
is held constant (fixed u, variable L). The figure indicates resonant behavior: for a given field strength, the binding 
energy passes through a maximum as L is increased from zero. 

The existence of a maximum binding energy follows from the limiting cases L <C R and L 3> R. The energy of 
the bound state is the solution of the secular equation, Eq. (|15p . When L is small, the binding energy increases in 
proportion to u 4 . When L is large, A(e) is proportional to L, and Eq. fT5"]) effectively becomes detGoTo = 0, which 
gives e = Aq. Corrections are of order RjL. Since the binding energy increases initially and approaches zero in the 
limit L 3> R, it must pass through a maximum. This leads to the peaks shown in Fig. [5l 

The decrease in binding energy with increasing length should not be interpreted as the uniform field limit of the 
square well model. The results of Fig. O were obtained by assuming the square well was smaller than the wavelength 
of a particle in the lowest band, which is not satisfied in a uniform field. 

These two figures indicate that a strong electric field applied over a small region of the nanotube is required to 
create localized states states far from the band edge. 
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O L.10 
L-20 
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Field Strength ( u ) 

FIG. 4: (Color Online) Binding energy as a function of potential strength in the square well model when the width of the 
potential is fixed. The horizontal axis shows the value of the dimensionless parameter u = eE^L / (Kvf / R) that characterizes 
the applied field strength. The lines are 7 = 1 — es/Ao for the potential widths L/R shown in the legend. The shape of the 
individual curves is similar to the binding energy curve of the delta function model in Fig. [3] Increasing the width of the well 
decreases maximum binding energy. 



V u = 5.0 
□ u = 2.0 
O U.1.0 
u = 0.5 




2 4 6 8 10 

Potential Width ( L/R ) 

FIG. 5: (Color Online) Binding energy as a function of potential width in the square well model when the field strength is 
fixed. The horizontal axis gives the width of the square well, L/R. The lines are 7=1 — e.g/Ao for the potential strengths u 
shown in the legend. The binding energy passes through a maximum as the width is increased, after which it falls back to zero. 



IV. CONCLUSION 



In this paper, we have considered the effects of nonuniform transverse electric fields on the electronic structure 
of semiconducting carbon nanotubes by using the linear Dirac Hamiltonian for states near the Fermi surface. We 
considered fields that vary slowly over distances comparable to the interatomic spacing, but rapidly over distances 
comparable to the wavelength of the envelope functions. This allows us to ignore intervalley scattering even though 
the effective potentials in the long-wavelength Hamiltonian change discontinuously. 

We studied two effective potentials: a delta function and a square well. In both models, the potential leads to the 
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formation of localized states inside the free nanotube band gap for arbitrarily weak fields. For weak fields, the energy 
of the subgap states moves away from the band edge in proportion to the fourth power of the electric field strength. 
When the range of the potential is fixed, the binding energy of the subgap states exhibits crossover behavior. The 
states lie close to the unperturbed band edge until the field strength exceeds a threshold value, after which they 
rapidly approach the center of the band. In the square well model, the binding energy of the subgap state passes 
through a maximum as the width of the well is increased. In both models, a strong electric field applied to a small 
region of the nanotube is required to create subgap states far from the band edge. 

The subgap states are exponentially localized around the center of the potential. As such a single subgap state is 
not likely to affect transport along the nanotube axis. However, these states might be important in processes involving 
transport perpendicular to the nanotube axis, such as scanning tunneling microscopy experiments. The effect on STM 
experiments designed to measure the nanotube band gap should be negligible, however. If the tip of the scanning 
probe just touches the nanotube, which just touches the substrate into which electrons will flow, then the electric field 
in the nanotube is roughly equal to the potential drop across its diameter divided by the static dielectric constant. In 
order for electrons to tunnel into the nanotube, the potential due to the STM tip must be on the order of the band 
gap: eAV » 2hv//3R. This gives u » 1/3. The subgap states at this field strength are virtually indistinguishable 
from the free nanotube band edge in either model. 

One might also consider several of these subgap states evenly spaced along the nanotube axis. This array could 
be used to control tunneling through the nanotube for charge carriers with energies inside the band gap. By tuning 
the energy of the subgap states, one could control tunneling across the nanotube from a reservoir of charge carriers 
at some energy e < Ao- Subgap states created by an external potential could also be important in exciton dynamics. 
They might provide specific sites along the tube for exciton localization and recombination. 

Experiments to explore this effect would require the generation of large electric fields — tens of volts per nanometer 
— in a relatively small region along the nanotube axis. One experimental signature of the subgap states would be an 
increase in tunneling current perpendicular to the nanotube axis with no change in transport along the tube axis at 
the same potential bias. 

The Green function approach we used is not limited to carbon nanotubes. The theory of graphene nanoribbons is 
similar to that of carbon nanotubes. The major difference is that periodic boundary conditions are replaced by hard 
wall boundary conditions. A uniform static electric field perpendicular to the edge of a nanoribbon applied over a 
small region would lead to a potential similar to Eq. ([7]). However, for nanoribbons, the applied field would mix all 
bands since the potential is linear instead of periodic. The methods described in this paper might be useful in spite 
of the differences. 

It would also be desirable to compare our predictions with the results of numerical simulations. However, it 
seems that the conditions in which our effective potentials accurately approximate the physical field make numerical 
simulation difficult. The potential breaks translation invariance and extends over a large number of unit cells, but 
a very small fraction of the length of the nanotube. As a result, the number of atoms and unit cells required would 
make a tight-binding simulation or density functional theory challenging. 

This work was supported by DOE grant No. DE-FG02-ER0145118, and by the NSF under grant DMR-00-79909. 
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APPENDIX A: COMPARISON OF ANALYTIC AND NUMERICAL RESULTS FOR UNIFORM STEF 



TableUcompares analytic and numerical predictions for the band gap response of a semiconducting carbon nanotube 
in a uniform transverse electric field. The first two columns give the chiral indices and radius of the nanotube analyzed. 
The third column gives the critical field strength predicted by Novikov and Levitovs^ 



E c w 0.62 • Hv F /eR 2 
205 



V/nm. 



The fourth column gives the approximate value of the threshold electric field strength from numerical simulations 
(the field strength at which the band gap started to decrease linearly with increasing field strength, interpolated from 
the graphs given in the references cited), and the fifth column gives the reference from which this value was taken. 
All numerical simulations were for semiconducting zig-zag nanotubes. The electric field strengths do not include 
depolarization effects. The values represent the field inside the nanotube, not the applied electric field. 



Tube 


Radius 


Critical Field 


Threshold 


Reference 


(m, n) 


(nm) 


(V/nm) 


(V/nm) 




(10,0) 


0.4 


2 


3.3 


Ref. .6 


(8,0) 


0.3 


3.2 


4 


Ref. W 


(11,0) 


0.4 


1.7 


2 


Ref. W 


(16,0) 


0.6 


0.8 


1.2 


Ref. 11 



TABLE I: Comparison of the critical field strength predicted by Novikov and Levitov and the threshold field reported in 
numerical tight binding calculations. The critical field from the linear Dirac theory provides a physical explanation for the 
crossover behavior observed in simulations. 

The Dirac Hamiltonian is the linear k • p approximation of a tight-binding Hamiltonian for carbon nanotubes, so 
the results of the two theories should agree, at least qualitatively. Novikov and Levitov acknowledge that higher- 
order terms break the symmetry that protects the band gapf£ so there is no true critical behavior. However, the 
critical field strength of the Dirac theory provides an intuitive explanation for the crossover behavior reported in 
numerical simulations. The low-energy Dirac theory captures the essential physics of carbon nanotubes necessary for 
a qualitative description of the response to an external electric field. 

APPENDIX B: NO DEPENDENCE ON CHIRAL ANGLE 

Here, we show that the density of states and the energy eigenvalues of the subgap bound states do not depend on 
the chiral angle. 

The Hamiltonian for a nanotube in a nonuniform transverse electric field can be written 

H(6) = e-^ e / 2 [(a y A m - iRa x d z ) S m , n + (u/2) f(z) S m , n±1 ] e^ ' 2 . (Bl) 
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This is equivalent to Eq. ([9]). (Since the potential is a scalar, e - lcr ^ 6 / 2 V m ^ n (z) e Icr * s / 2 = V m . n {z).) We will refer to 
the Green functions derived from this Hamiltonian as "chiral." Those studied in the main body of the text do not 
depend on the chiral angle, so we will refer to them as "achiral." 
The equation satisfied by the chiral Green functions is 

[eS m ,n ~ Hm,n(«)] ©m(0) = (B2) 

Since e is also a scalar, we may write 

e-**'" 2 (sS m , n - H m , n ) e™* 6 ' 2 6 m (6) = S(z). (B3) 

Multiplying on the left by e la * s l 2 and on the right by e - %a ^ e / 2 ^ we obtain an achiral operator acting on a transformed 
Green function. The transformed Green function is the inverse of an operator that does not depend on the chiral 
angle, and so it too must be independent of chiral angle. The equation satisfied by the achiral Green functions is: 

(e<W - H m , n ) & m = S{z). (B4) 

The chiral and achiral Green functions are related by a unitary transformation: 

& rn = e i °* 9 l 2 & rn {6)e- i °* e ' 2 . (B5) 

The results of the main body of the paper apply to the achiral Green functions. Next, we show that the bound 
state energies and density of states derived from achiral Green function & m are the same as those derived from the 
chiral Green functions. This implies that the density of states and bound state energies do not depend on the chiral 
angle. 

The density of states depends only on the trace of the Green function. Since the trace of a matrix is invariant under 
unitary transformations, the density of states derived from 6 m and & m {9) are identical. 

The energies of the bound states of the m-th band are the poles of the chiral Green function, which occur when the 
determinant of the denominator in Eq. (|13|) vanishes: 

det{l - G m T m } = 0. (B6) 

The matrix r m is a linear combination of two achiral Green functions: 

T m = aG m+ i + /?G m _i. (B7) 

As a result, the relation between T m and the chiral matrix T rn (9) is identical to that between the Green functions, 
given in Eq. (IB5|) . 

Since the identity matrix is invariant under unitary transformations, Eq. (|B6j) can be expressed in terms of the 
chiral Green functions as 

det {e"^ 2 [1 - G m (e)r m (6)] e^/ 2 } = 0. (B8) 

The determinant of a matrix is invariant under unitary transformations. Therefore, the determinants of the chiral 
the achiral Hamiltonian are equal, and lead to the same electronic spectra. 

To summarize, neither the density of states nor the energies of the subgap states depends on the chiral angle of 
the nanotube for the models studied in this paper. This is consistent with known results for a free nanotube or a 
nanotube in a uniform transverse electric field. In these models, there is no dependence on the chiral angle when the 
effective Hamiltonian includes only terms linear in d z and d^, (or q and A m in Fourier space). 

There are two ways in which a dependence on the chiral angle might be introduced. First, one could introduce 
the quadratic terms of the k • p approximation to the effective Hamiltonian with a scalar potential. Trigonal warping 
cannot be removed with a unitary transformation and leads to a dispersion relation that depends on the chiral angle. A 
second possibility is to introduce a potential that mixes states on different sublattices. Even if the effective Hamiltonian 
is linear in d z and d^, the potential will introduce a term proportional to a x , a y , or some linear combination of the 
two. In this case, the chiral Hamiltonian TL{9) is no longer related to an achiral operator by a unitary transformation, 
and the eigenvalue spectrum will depend on the chiral angle. 
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APPENDIX C: DERIVATION OF GREEN FUNCTIONS 



Here, we derive the achiral Green functions for the unperturbed Hamiltonian. As discussed in App. [B] the bound 
state energies and the density of states do not depend on the chiral angle, so we ignore it to simplify the calculation. 

Because the free nanotube Hamiltonian is static and invariant under translations along the nanotube axis, the 
Green functions can only depend on the difference between coordinates: G(z, z') = G m (z — z'). We set the z' to zero 
when solving for the Green function. The equation satisfied by G m (z) is 

(e - (T y A. m + iRa x d z ) G m (s, z) = S(z). (CI) 

This linear differential equation can be solved by Fourier transformation: 

G m (e, z)= [°° p- e-*' G m (e, q). (C2) 

Using this expression for G m (e, z) in Eq. (|C1|) . the solution is 

^ , s £ + cr y A m + cr x qR 

Gm(£ ' g) = ^-( g fl)'-(A m)a + ^ - (C3) 

A convergence factor ir] has been added to the denominator to push the poles of the Green function into the appropriate 
quadrants of the complex plane. This Green function will satisfy boundary conditions for outgoing waves as z 
approaches ±00. 

Since the potentials considered in this paper break translation symmetry along the nanotube axis, the real space 
Green function is easier to use in calculations. Inverting the Fourier transform for q gives 

G m (e,2)= p-e*"G m {e,q) 

J-00 27T 

_ s + <r y A m + sgn(z - z 1 ) a x C m (e) iU(e) \ Z „ Z >\ /R (r . 



where £ m ( £ ) = \J s 2 — A^. This is equivalent to the Green function given in Eq. (fl 



APPENDIX D: JUSTIFICATION FOR IGNORING TRANSITIONS TO BANDS WITH Am| > 1 

Here, we argue that the error introduced by ignoring all terms in the Dyson series that involve transitions to bands 
other than the two nearest neighbors is small for energies near the band edge. 

To keep the expressions simple, we consider the full Green function for z — z' = and suppress dependence on z, 
z' ', and e. In the delta function model, the series for the full Green function in the lowest band is 

S = G + u 2 G (G+i + G_i) G 

+ m 4 Go(G+i+G_i) Go(G+i+G_i) Go (Dl) 
+ m 4 Go (G + iG +2 G+i + G_iG_ 2 G_i) G + 0(u 6 ). 

All terms of the infinite series involving only transitions to the two neighboring bands can be generated from a 
single term, which we call S: 

E = u 2 (G+i +G_i) G . (D2) 

The 0(u 2n ) term of the series can be expressed as (S) n plus other terms that involve Green functions from bands 
other than m = 0, ±1. For instance, the 0(u 4 ) terms can be expressed 

G S 2 + u 4 G (G+iG +2 G +1 + G_iG_ 2 G_i) G . (D3) 

In the delta function model, the terms of the series are products of the free nanotube Green functions in Eq. (I10p 
evaluated at z = z' = 0: 



14 



Near the band edge, e « Aq. Therefore, the denominator of Go is very small, but the denominator of G n with ti^O 
is of order unity. As a result, the largest terms of the series at any order are those with the most factors of Go in the 
product — i.e. the terms generated by (£)". 

As an example, we compare the two 0(u ) terms 

GoG+iGoG+iGo ~ | ; 1 | [ 1 | (D5) 

and 

We have ignored all common factors, as well as the numerators, which are matrix products of terms of order unity. 
Eq. (|D5|) is generated by (X) 2 . It is greater than the second term by a factor of y/e 2 — (A+2) 2 / \J e 2 — (Ao) 2 , which is 
large when e is near the band edge. For instance, if |e|/|Ao| = 0.1, the first term is large than the second by a factor 
of roughly 16. 

Thus, the error introduced by discarding all terms except those generated by S is small when the energy is near 
the band edge. We approximate the full Green function by a partial summation 

oo 

6o«G ^(S)". (D7) 

n=0 

Multiplying by the right by E, we have 

oo 

6 E-G ^(E) n+1 = 6o-G . (D8) 

n=0 

This can be solved to give 

6 = Go(l-£r 1 . (D9) 

(Formally, this series only converges when the eigenvalues of £ are less than onei 17 i 18 ) 

Eq. (|14[) gives the most general solution, with z ^ z' , and its derivation closely parallels the steps in this appendix. 
If z = z', Eq. (pi]) reduces to G (1 - E)" 1 . 

APPENDIX E: BOUND STATE FOR ANY FINITE POTENTIAL 

Here, we show that subgap states exist for arbitrarily weak fields in the delta function model of Sec. |TTJ 
If there are poles inside the band gap, then the denominator of Eq. (|14p must vanish for some energy in the range 
< e < Ao. To simplify the expressions that follow, we define C = e/Ao as the ratio of the particle energy to the 
unperturbed band gap, and x — uL/AR as a dimensionless measure of the potential strength. Poles will be solutions 
of the equation 



= det 



a _ x 2 ( + °y I C + Ml + l/Ao) , C + a y (l-l/A ) 



VW^ W(l + l/A ) 2 -C 2 J(l-1/A Q ) 2 -C 2 , 



(El) 



After multiplying out all the terms, the resulting product is a matrix of the form M = A(() 1 + B(()a y . The 
determinant of such a matrix is A(£) 2 — B(Q 2 . In this case, substituting Aq = 1/3 gives 



2 (C 2 + 4) v^e + (C 2 - 2) yie^c 2 



A(C) = 1 - x 2 " " ' , (E2) 

V(l-C 2 )(16-C 2 )(4-C 2 ) 



, 2 ^16^-5^4^ 



B(C) = x , V (■ (E3) 

V(W 2 )(16-C 2 )(4-C 2 ) 



15 



If there are bound states inside the band gap, then we expect their energy to be close to the band edge when the 
potential is small i.e. if uL <C R, then £ ss 1. Letting ( = 1 — r\ then Taylor expanding A(() 2 — B(() 2 = to 0{rj) 
gives 

= (90 - 36\/5 ^3 - Vtj X 4 ) V - 30V6 (VE - lj X 2 Vv + 0{rf' 2 ). (E4) 

Neglecting the terms of order if' 1 2 and higher, this equation has two solutions. One of these is r\ = 0. This solution 
corresponds to e = Ao, which is not a pole of the full Green function in Eq. ([H|) . The other solution is 

V = (^-i) (E5) 

where r ~ 0.991 and s w 0.677. Thus, any nonzero applied field produces a bound state inside the band gap, and 
the binding energy of this state scales with the fourth power of the potential strength when the potential is small. 
The singularity in this expression is a consequence of truncating the Taylor series. The expression diverges when 
X ~ 1.01, but the expression was derived assuming x "C 1. Fig. [3] shows that the energy of the subgap states has no 
singularities. 



